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' Abstract 

' Wc investigate mass losses via stellar winds from sun-like main sequence stars with a wide range of 

, activity levels. We perform forward-type magnetohydrodynamical numerical experiments for Alfvcn wave- 

fT^ ' driven stellar winds with a wide range of the input Poynting flux from the photosphere. Increasing the 

magnetic field strength and the turbulent velocity at the stellar photosphere from the current solar level, 
the mass loss rate rapidly increases at first owing to the suppression of the reflection of the Alfven waves. 
The surface materials are lifted up by the magnetic pressure associated with the Alfven waves, and the 
cool dense chromosphere is intermittently extended to 10 ~ 20 % of the stellar radius. The densities of 
' the corona and transition region above the chromosphere is also high, which leads to efficient radiative 

Qh, losses. Eventually most of the input Poynting energy from the stellar surface escapes by the radiation. As 

Q ' a result, there is no more sufficient energy remained for the kinetic energy of the wind; the stellar wind 

( , saturates in very active stars, as observed in Wood et al. The saturation level is positively correlated with 

c/5 ' Br.ofo, where Br,o a-nd /o are the magnetic field strength and the filling factor of open flux tubes at the 

^ C^^l photosphere. If /o is very small, namely a stellar surface is mostly covered by closed magnetic loops, and 

Br,ofo is small ^ 0.5 G, the possible maximum of the wind kinetic energy is comparable to that of the 
current sun. On the other hand, if i3r,o/o is relatively large ^ 5 G, the mass loss rate could be as high as 
J> ■ 1000 times. If such a strong mass loss lasts for ^ 1 billion years, the stellar mass itself is affected, which 

I could be a solution to the faint young sun paradox. Our simulations show that 0.1 - 10 % of the input 

Poynting flux is transferred to the final kinetic energy of the stellar winds. We derive Reimers-type scaling 
relations that estimate the mass loss rate and the wind kinetic energy of sun-like stars from this energetics 
consideration. 

■ Key words: magnetic fields - stars: coronae - stars: late type - stars: mass loss - stars: winds, 
'. outflows - waves 



_ ^ 1. Introduction temperature is also higher (Ribas et al. 2005; Telleschi et 

' al. 2005). A young stars have a very thick chromosphere 

. Main sequence stars with mass roughly below 1.5 so- extending to 10-20% of the stellar radius (Czesla et al. 

lar mass {Mq hereafter) posses a surface convective layer, 2012), which is in contrast to the thin chromosphere of 

which is the main origin of X-ray and wind activities the present Sun with width of 0.1-1% of the solar ra- 

from such solar-type stars; magnetic fields are gener- dius (Rq hereafter). These observations show that in 

ated in a surface convective layer by dynamo mechanisms young sun-like stars the atmospheric materials are lifted 

(e.g., Choudhuri et al.l995; Brun et al.2004; Hotta et up to higher altitudes and heated up more intensely in 

al.2012) and turbulent motions associated with the mag- upper regions. This is probably related to the strong 

netoconvection are the source to drive flares and outflows, magnetic fields with an order of kG or even larger ob- 

Such magnetic activities are themselves interesting phe- served in young main sequence stars (Donati & Collier 

nomena with fruitful physics, and in addition they may Cameron 1997; Saar & Brandenburg 1999; Saar 2001; see 

also affect planetary circumstances around central stars, also Donati & Landstreet 2009 for recent review), which 

Various researches have been carried out recently from this are much stronger than the average strength of 1-10 G of 

viewpoint (e.g., Terada et al.2009; Stercnborg et al.2011; the present-day Sun. 

Lammer et al.2012). On the other hand, the observed mass loss rates by 
Young sun-like stars are very active: The observed X- stellar winds show a different trend. Observations of near- 
ray flux is up to ^ 1000 times larger than the present so- by stars with the spectral types of G,K,M by Wood et 
lar level (Giidel et al. 1997; Giidel 2004), and the X-ray al.(2002; 2005) show that, after an increase of the mass 
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Fig. 1. Geometry of a flux tube and the input parameters. 

loss rate with the increasing X-ray flux, it saturates at ^ 
100 times of the current solar level and even drops in some 
very active stars. An explanation of the saturation is due 
to the change of the magnetic topology. It is expected that 
the surface area of active stars is mostly covered by closed 
magnetic structure as a result of the strong magnetic fields 
(Schrijver & Aschwanden 2002). Then, the atmospheric 
gas is confined in closed loops rather than streaming out 
from open flux tubes (Wood ct al. 2005; Vidotto et al. 
2009). 

In this paper, we propose an alternative and additional 
mechanism for the saturation. We focus on the dynamics 
in open magnetic flux tubes and show by numerical simu- 
lations that the saturation of the stellar winds could take 
place in open flux tube regions because of an enhancement 
of the radiative losses in active stars. The construction 
of the paper is as follows: After briefly introducing the 
setup of the MHD simulations in §2, we summarize in §3 
basic formulations for the analysis of the wind energetics 
in the later sections. §4 presents the main results of the 
simulated stellar winds. In §5 we extensively discuss our 
results in light of observations of stellar winds (§5.1) and 
chromospheres (§5.3). In particular, we directly compare 
the simulated winds with the observations by Wood et 
al.(2002; 2005) in §5.1. In §5.2, we discuss the resuhs in 
terms of the faint young Sun paradox. 

2. Simulation Setup 

2.1. Setup 

We extend our one dimensional (ID) magnetohydrody- 
namical (MHD) simulation code originally developed for 
the present-day solar wind in Suzuki & Inustuka (2005; 
2006) to young active sun-like stars with strong magnetic 
fields and large velocity fiuctuations at the photosphere 
(§2.2). In this paper, we focus on main sequence stars and 
fix the basic stellar parameters on those of the present-day 
Sun, mass, M = Mq, radius, rg = Rq and effective tem- 
perature. Toff = 5780 K. Compared to the current Sun, the 
radius of a IMq star is slightly smaller by ^ 10 % and the 



effective temperature also slightly lower within < 200 K at 
early epochs of the main sequence (Sackmann et al. 1993). 
However, we use the same radius and effective tempera- 
ture in our simulations for both active and inactive stars 
because we would like to study responses of the stellar 
winds to the changes of the magnetic fields and velocity 
fluctuations at the photosphere^. 

We dynamically solve ideal MHD equations with ra- 
diative cooling and thermal conduction in super-radially 
expanding flux tubes from the photosphere at r = rg to the 
outer boundary at r ss 307'o. We only consider the deriva- 
tive with respect to r but treat the three components of 
velocity and magnetic field. We initially set up a static 
and cool (T = 10** K) atmosphere and start simulations by 
injecting the three components velocity fiuctuations from 
the photosphere with amplitude of each component, Svq, 

{5vl)^ T'^^ P{u:)du:, (1) 

where we adopt a frequency spectrum, P{uj) oc from 
l/i^max = 20 seconds to 1/wmin = 30 minutcs. The nor- 
malization of P{uj) is determined to give Svq. 

Various modes of MHD waves are excited by the surface 
fluctuations (e.g. Isobe et al.2008; Kato et al.2011), and 
propagate upwardly (e.g., Bogdan et al.2003; Okamoto & 
De Ponticu 2011). Propagation, reflection, and dissipa- 
tion of these MHD waves are automatically treated. The 
Alfven wave, among various modes, is supposed to play a 
major role in driving stellar winds mainly because it trav- 
els a long distance to the wind acceleration region (e.g., 
Alazraki & Couturier 1971; Hollweg 1973). A great advan- 
tage of the code is that we can determine mass loss rates 
as a direct output of the injected Foynting flux from the 
photosphere. On the other hand, we cannot handle waves 
with wavelengths shorter than the grid scale. These waves 
dissipate and heat up surrounding gas to conserve the to- 
tal energy equation in our treatment. In other words, 
we implicitly assume that sub-grid-scale structures are in- 
stantly transferred to the thermal energy by cascading 
processes or shocks. 

As for the cooling in the coronal region, we adopt the 
cooling table for the optically thin plasma with the solar 
abundance (Landini & Monsignori-Fossi 1990; Sutherland 
& Dopita 1993), in which the cooling rate, qr, (erg 
cm^'^s^^), is proportional to . For the cooling in 
the chromosphere, we adopt an empirical cooling rate 
gj^ = 4.5 X lO^p erg cm~^s~^ by Anderson & Athay (1989). 
This cooling rate takes into account an effect of opti- 
cally thick cooling under the non-LTE conditions (but 
see Carlsson & Leenaarts 2012 for a more detailed treat- 
ment based on snapshots of 2D MHD simulations). We 
smoothly connect the cooling rates in these two regimes 
by an interpolation. In the transition region, the cooling 
rate is initially cx p at lower heights but shifted to at 
higher altitudes. 

^ Also, sucii subtle changes of the radius and effective temperature 
give almost negligible effects on the dynamics and the energetics 
of the stellar winds. 
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Table 1. Labels for input parameters. For example, the label 
for the standard case is 'CbO«'. 



We do not explicitly take into account stellar rotation in 
the dynamics of the stellar winds. If a star rotates by more 
than ~ 20 times faster than the present Sun, the termi- 
nal velocity of the wind will be affected (§5.4). However, 
the mass loss rate is affected little (Belcher & MacGregor 
1976). The stellar rotation rather significantly affects the 
strengths of the generated magnetic fields (Holzwarth & 
Jardine 2007). In this paper, we consider this effect by 
incorporating wide ranges of the parameters on the mag- 
netic fields. 

2.2. Input Parameters 

We investigate how the properties of the stellar winds 
depend on the following four parameters (Figure 1): 
Velocity perturbation, Svq., and radial magnetic field 
strength, Brfi, a filhng factor, /o, of open^ flux tubes, 
which are all measured at the photosphere, and a typical 
height, hi, of closed loops, which surround an open flux 
tube we are considering. h\ determines the location of the 
rapid expansion of a flux tube (Figure 1). From /o and 
hi, we adopt a functional form of a filling factor'^ which 
depends on r. 



.fir) 



e^^+/o-(l-/o)/e 



(2) 



1 



This is the same functional form as in Kopp & Holzer 
(1976) with adopting h\ = Ri — Rq = a in Equation (11) 
of their paper. The conservation of magnetic flux fixes 
radial magnetic field. 



Br — Br 



f(r)r' 



(3) 



As a reference case for the present-day solar wind, we 
use the following values, which explain observed proper- 



We call an 'open' flux tube if the magnetic fields do not close in 
the simulation region. 

In many literatures including ours (Suzuki & Inutsuka 2005; 
2006) a super-radial expansion factor, instead of a filling factor, 
is used to set up open flux tubes. Defining /'(r) as a super- 
radial expansion factor, these two factor are simply related as 

/(0 = /o/'M. 



ties of polar coronal holes: (5fo = 1.34 km s ^, Br,o = 1 
kG, /o = 1/800, and h = O.Oli?0(= 7 x 10^ km).' The 
value of 5vq is within the range of observed granulations 
(Holweger et al. 1978; Matsumoto & Kitai 2010). We 
adopt the values of B^.o and /o to explain recent ob- 
servations by HINODE (Tsuneta et al. 2008; Shimojo 
& Tsuneta 2009; Ito et al. 2010; Shiota et al. 2012). 
In these papers, the authors reported that there are a 
number of super-radially open flux tubes anchored from 
strong magnetic field patches with B^.o ^ kG at the pho- 
tosphere in polar regions. Comparing the photospheric 
magnetic field strength (^ 1 kG) and measured interplan- 
etary field strengths normalized at the earth orbit (^1 — 10 
uT = 10-5-10-4 G; e.g.. Smith & Balogh 2008), a typical 
filling factor, /o, of open field regions can be estimated as 
an order of 1/1000 (see Equation 3). 

We perform 163 runs in wide ranges of the parameters: 
5vq = (0.669 - 7.57) km s-^, B^.o = (0.5 - 16) kG, /o 
(1/400-1/6400), and hi = (0.01 - 0.1)Vo. In active young 
stars, Brfi is stronger, h\ is larger (taller loops), while 
closed loops are more dominant than open field regions, 
and therefore, /o is smaller. We naively expect that 5vq is 
larger in active stars, whereas too strong magnetic fields 
might inhibit footpoint motions of flux tubes (Katsukawa 
& Tsuneta 2005). We label the simulated models by (A 
- F) for Brfifo, (a - f) for i?r,0: (—2 — 1-5) for Svq, and 
(a - 7) for h\, which are tabulated* in Table 1. Using 
these labels, for instance the standard case (Sr.o./o = 1-25 
G, Br,o = 1 kG, dvo = 1.34 km s-\ and hi = O.Olro) is 
labeled as 'CbOa'. 

3. Energetics Formulation 

In this paper, we investigate how the mass losses are 
controlled by the surface properties by examining the en- 
ergetics of the simulated stellar winds. We flrstly consider 
the time-averaged structure of each simulation run in the 
next section (§4). In this section we summarize some ba- 
sic equations for the wind energetics. Under the steady 
state conditions, an equation of the total energy can be 
written as (e.g., Fisk et al.l999; Suzuki 2006) 

7 p GM 



7- 1 p 



47r 



(tj X B) X J5 + Fc 



9R = 0, 



(4) 



where is thermal conductive flux, and the other vari- 
ables have the conventional meanings. We focus on the 
radial component of the equation from now. We can write 
the divergence of an arbitrary vector. A, as 

1 d 



'-f dr 



[r'fAr 



(5) 



to take into account the super-radial expansion of an open 



The labels for Svq are the twice of the power-law indices of 2 
scaled by Svq = 1.34 km of the standard case. For example, 
'-I-3' for Svo = 3.79 km s'^ is from (3.79... )2 = 2+^ x (1.34.. .)2 
or 3.79... = 2+3/2 1.34^^^. 
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flux tube. The integration of Equation (4) from an arbi- 
trary reference point, rrof, to r gives 



Airr f 



= Ml- 



7- 1 p 
, 1 P . 



7 p GM\ 



■An 



47r 



7 — 1 p A-Kp 



-47rr^/Fc,r + 47r / gRr^/d?' = const., (6) 



where subscript, _L, indicates the perpendicular compo- 
nents to the radial direction, M is mass loss rate, 



i'l = Anr'^fpvr, 
and $B is radial magnetic flux, 
$B=4^r2/Br. 



(7) 



(8) 



The energy flux of Alfven waves along with r direction 
can be written as 



^^-^^l^f +4^ -^'T 



(9) 



(Jacques 1977; Cranmer ct al. 2007). One may notice 
that Fa partly consists of Equation (6). Integrating over 
the total area of the open flux tubes on a stellar surface, 
we can define Alfven wave luminosity. La, as 



LAir)f{r) 



47rr^/(r)LA(r 
2 S2 



(10) 



Airp 



v±Bj 
47r 



Since we consider the energy transfer in super-radially 
open flux tubes, luminosity, L^f, from open flux tubes, 
instead of La integrated over the total area, should be 
treated in the energy conservation. 

The right-hand side of Equation (9) is separated into 
the two parts; the first term indicates the energy flux ad- 
vected by background flow, Vr, and the second term, which 
exists even in static media, is the Poynting flux involving 
magnetic tension. Wc can introduce Elsasser variables, 

z± = v±T^=, (11) 



where z+(z_) denotes an amplitude of Alfven waves which 
travels to the (anti-)parallel direction with Br. Using z±, 
the second term of Equation (9) can be expressed as 



— Br ~ 



Si 



1 



4tt 



p(z\~zi)vA, 



(12) 



where 7JA - 



is Alfven velocity. This expression illus- 
trates that \pz\ [jpz"^ ) corresponds to the energy den- 
sity of Alfven waves to the (anti-)parallel direction with 
Br, and that —Br ^^^^ corresponds to the net positively 
traveling (parallel with Br) Poynting flux associated with 
Alfven waves in static media. 



After the injection of Alfven waves from the photo- 
sphere La/ decreases because of various types of damp- 
ing processes, such as turbulent cascade (Matthaus ct 
al. 1999; Chandran 2005; Cranmer et al. 2007; Verdini 
& Velli 2007), phase mixing (Heyvaerts & Priest 1983; 
Grappin et al. 2000; De Moortel et al. 2000), nonhn- 
ear mode conversion to compressive waves by magnetic 
pressure, B\/%n (ponderomotive force; Kudoh & Shibata 
1999; Suzuki & Inutsuka 2005; 2006) and by paramet- 
ric decay (Goldstein 1978; Terasawa et al. 1986; Suzuki 
& Inutsuka 2006; Nariyuki 2009), ioncyclotron reso- 
nance (e.g., Axford & McKcnzic 1997; Tu & Marsch 
2001). Among these processes, our ID MHD simulations 
can only consider the nonlinear mode conversion to com- 
pressive waves. Recently, Matsumoto & Suzuki (2012) 
have performed 2D MHD simulations, which also handle 
phase mixing and a part of turbulent cascade. This result 
gives similar wind structure to the ID result by Suzuki & 
Inutsuka (2005). Therefore, we think that the ID simu- 
lation, which we use in this paper, still gives reasonable 
properties of the stellar winds, even though the actual dis- 
sipation channels might be different. After the damping 
of the Alfven waves, a small fraction of the input energy 
is flnally transferred to the kinetic energy of stellar winds. 



LK(r)/(r)-M^=47rrV(r)p«.|-^-' 



47rrV(0-F^K(r),(13) 



where Lk and Lk are kinetic energy flux and luminosity. 

Our main aim is to understand how the wind ki- 
netic energy luminosity, LK(rout)/(?'out) =LK(rout) (note 
fifout) = 1), at the outer boundary, rout ~ 30ro, is deter- 
mined by the input Poynting flux (La/)o at the photo- 
sphere. Hereafter, we express, e.g. (La/)o and LK,out, 
etc. for LA('ro)/o and LK(?'out) for simplicity. In our 
simulations, we only inject velocity perturbations with 
6vo from the photosphere without magnetic field pertur- 
bation. Then, the input energy flux can be written as 
La,o = Po{^Vq)va,o in the static background {vr =0), where 
(Svq) is kinetic energy per mass of the sum of the two 
transverse components. Then, the input wave luminosity 



(La/)o = 47rro/opo(<5wo)^'A,o- 



(14) 



Wc examine the energetics of the stellar winds by 
two steps (Figure 2); we firstly inspect the energetics 
in the chromospheres. Here, the reflection of outgo- 
ing Alfven waves play an important role because the 
Alfven speeds change rapidly because of the steep decrease 
of the densities. We define r = rtc as the top of the chro- 
mosphere where T = 2 x 10^ K, and wc measure (LA/)tc 
at r = rtc to quantify the reflection (§4.2). We inspect the 
transmissivity. 



Ct = {LA.f)tc/{LAf)o, 



(15) 



of the Alfven waves through the chromosphere in different 
cases. 

Next, we examine how much fraction of the surviving 
wave energy, (LA/)tc, is finally transferred to the wind 
kinetic energy (§4.3). By comparing dominant terms at 
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Fig. 2. Schematic picture for the energetics in an open flux 
tube. Please see text for the notation of each variable. 



Ttc and at Tout in Equation (6), wc can derive an energy 
conservation relation: 



where 



(-^R/)tc — {LcDtc, 



{Lnf)t 



r°"' 2 

47r / gRT fdr 



(16) 



(17) 



is the radiation loss from r ~ rtc to the outer region, and 
the gravitational loss, 

[L^Dtc^M—. (18) 

Please note that the energy loss by the downward ther- 
mal conduction in the transition regions is included in the 
radiation loss term, (LR/)tc, because the thermal con- 
ductive flux from the upper coronae mainly escapes by 
the radiation in the transition regions. 



4. Results 

Photosphere 



Wind connection 



Before carrying out detailed analyses, we show the rela- 
tion between the input Alfven wave luminosity {LAf)o 
from the photosphere and the output kinetic energy, 
iK,out. The upper panel of Figure 3 shows the (La/)o ~ 
-^K,out diagram. In the figure we overplot three lines 
for constant conversion factors, ce = 0.001,0.01,0.1, from 
(La/)o to LK,out by the three lines. Using Equations (3) 
and (14) with va,o = Br,o/V^''^PO: we can explicitly write 
the relation that connects iK,out to {Lxf)o'- 



^K.out = C^{Lp,f\) = CE$B\/ — ((^Wi 
f Br,ofo \/f Svq 

V 1-25G / \ I 1.34km s-i 




10^'^g cm^-^ 



1/2 



(19) 
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Fig. 3. top: Relation between the final kinetic energy lumi- 
nosity, Z/K.outi a-nd the input wave energy luminosity, (La/)Oi 
at the photosphere of the simulated stellar winds. The lines 
indicate constant conversion factors, cb = 0.1,0.01,0.001. The 
right and top axises are normalized by the solar luminosity, 
Lq = 4 X lO'^'^ erg s~^. bottom: Relation between M and 



(La/)o divided by 
conversion factors, cm : 



- = iiM The lines indicate constant 
: 0.1,0.01,0.001. 



which is normalized by the input parameters and the out- 
put LK,oout = 2.1 X 10^^ erg with Ce =0.017 of the stan- 
dard case for the present-day solar wind (Model CbOa). 

Different colors are used for different sets of Sr,o/o- We 
would like to note that i3r,o/o determines the magnetic 
field strength in the outer region where the super-radial 
expansion of the flux tubes is already completed { f{r) — >■ 1 
in Equation 3). Sr.o/o is a good indicator of the properties 
of the solar wind, particularly the wind speed (Kojima 
et al. 2005; Suzuki 2006). The figure shows that the 
conversion factor is distributed in a range of 0.001 < ce < 
0.1, and a typical value is ce ^ 0.01. 

Although the scatters are not small, focusing on a sin- 
gle set of Brfifo, Ce shows a trend; ce initially increases, 
namely LK.out increases with (La/)o much faster than the 
linear, which is followed by a decrease of ce or a satura- 
tion (or even decrease) of iK,out in the large (La,/)q- The 
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initial increase of ce can be explained by the suppression 
of the reflection of the Alfven waves (Suzuki 2012), which 
we discuss in §4.2. The later decrease is due to the en- 
hanced radiative losses (§4.3). The figure shows that the 
saturation level of i^K.out is determined by i?r,o./o: which 
we also discuss in §4.3. 

It is more useful for readers to show the mass loss rates 
as a function of surface properties. The lower panel of 

Figure 3 shows M with wave luminosity divided by = 
where Wcsc,o(= 618 km s^^) is the escape velocity. 
Both axiscs arc plotted in unit of A/0yr~^. Similarly to 
the scaling relation for LK.out (Equation 19), we can define 
a conversion factor, cm, with respect to the mass loss rates. 
Then, we can explicitly write down an equation that gives 
M of general sun-like stars with radius, tq, and mass, M, 
from the surface properties in a similar manner to the 
Reimers (1975) formula (sec also Schroder & Cuntz 2005) 



M- 



(iA./)o^o 




20) 



where the normalizations are adopted from the input pa- 
rameters and the results {M = 2.2 x lO"^^M0yr~^ and 
Cm = 0.023) of the standard case for the present-day Sun. 

In general, density, po, at a photosphere can be esti- 
mated from the stellar basic parameters. For instance, in 
the cool star condition with effective temperature, 4000 

K< Teff < 6000 K, Po cx {j^^ T^a (§9 of Gray 1992; 

Suzuki 2007). Equation (20) is a general formula to es- 
timate mass loss rates of stars with a surface convection 
zone (see also Cranmer & Saar 2011). Once the basic stel- 
lar parameters, rg, M, and T^s, are obtained, Equation 
(20) has the three undetermined parameters, cm, -Br,o/o, 
and Svq. One can use a typical value, cm ~ 0.02 from 
the simulation results. Svq is expected to be a fraction of 
the sound speed at the surface; Svq = 1.34 km s~^ of the 
standard case corresponds to '--^ 20 % of the sound speed. 
Br.ofo, magnetic field strength in open field regions, is the 
most unknow parameter among the three. If M is ob- 
tained by other methods. Equation (20) can be used to 
derive i?r, o/o in the opposite way. 

The lower panel of Figure 3 shows that the average 
trend of cm is slightly increasing with while the 

upper panel illustrates that the average trend of ce does 
not show a systematic trend with (La/)o- This is mainly 
because the wind speed becomes slower as the density in- 
creases with increasing (La/)o; the dense wind cannot be 
accelerated to high speed. If a star rotates very rapidly, 
the wind speed of active cases (on the right side of the 
panels) might be higher. In this case ce might show an 
increasing tread with (iA/)o, similarly to cm. 
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Fig. 4. Surviving fraction, ct = (iA/)tc/(-f'A/)o, of the 
Alfven wave luminosity at tlio top of the chromosphere versus 
input wave luminosity at the photosphere, (La/)o. 



4.2. Initial Increase of Lk,i 



-Wave Reflection- 



4-. 2.1. Energetics 

A sizable fraction of the input Alfven waves suffers re- 
flection before reaching the corona (Hollweg 1984; Moore 
et al. 1991; Murawski & Musielak 2010; Verdini, Grappin, 
& Velh 2012). Obsevation by HINODE obtained sig- 
natures of reflected Alfven waves at the photosphere 
(Fujimura & Tsuncta 2009). In the chromosphere, the 
density decreases very rapidly, which leads to the rapid 
increase of the Alfven velocity. Therefore, the shape of 
outwardly propagating Alfven waves is considerably de- 
formed, which indicates that these waves are reflected 
back downward. For the current solar condition, typically 
^ 90 % of the input energy is reflected back downward 
(Suzuki & Inutsuka 2006; Matsumoto & Suzuki 2012). 

We measure (iA/)tc at the top of the chromospheres, 
r = rtc, to quantitatively study the surviving fractions 
(= Ct in Equation 15) after the propagation through the 
chromopshres. The reflected waves are mostly dissipate 
through the downward propagation in the chromospheres 
and finally lost by the radiation, where some fractions 
are considered to be again reflected upward and may be 
trapped (Matsumoto & Shibata 2010). Thercforcm the 
radiation loss in the chromsphcres is mostly taken into 
account in the component of the reflected waves. 

Figure 4 displays ct = (iA/)tc/(iA/)o- For small 
(La/)o, the reflection of the Alfven waves is so effective 
that only ^ 1% of the input energies can transit through 
the chromospheres. On the other hand, the fraction is in- 
creasing with increasing (La/)o and finally reaches ^ 0.3. 
This trend can be understood from the density structures 
in the chromospheres, which is discussed later (Fig. 5). In 
short, in cases with large energy inputs the reflection of 
the outgoing Alfven waves is not so effective because the 
densities decrease more slowly and accordingly the change 
of the Alfven speeds is more gradual. 
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Model 


SVQ 


Br,o(kG)x/o 


hi (ro) 


(LaDo 


(iA/)tc 


CT 


M (Moyr-i) 


f r.out 




(iR/)tc 


Ba-2a 


0.669 


0.5/800 


0.01 


1.5 X 10^» 


5.8 X 10^" 


0.038 


1.3 X lO"-''"' 


777 


2.4 X 10^" 


1.2 X 10^" 


CbOa 


1.34 


1/800 


0.01 


1.2 X 10^'' 


6.4 X 10^' 


0.052 


2.2 X 10"-"* 


554 


2.1 X 10^' 


1.5 X 10^ ' 


Cd+37 


3.79 


4/3200 


0.1 


9.8 X 10^^ 


1.6 X 10^" 


0.17 


8.5 X 10--*-' 


313 


2.6 X 10^** 


2.2 X 10^» 



Table 2. Input parameters and outputs for the cases shown in Figure 5. The unit of Svq and Vr^out is km s ^, and the unit of L is 
erg s^^. ct is the tranmission fraction of the Alfvcn waves from the photosphere to the top of the chromosphere (Equation 15). 
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Fig. 5. Comparison of the time-averaged wind structures of 
the three cases tabulated in Table 2. Dashed, dotted, and 
solid lines are the results of Models Ba— 2a, CbOa (the stan- 
dard case for the present-day Sun), and Cd+S-y. From top to 
bottom, the densities, radial velocities, and temepratures are 
compared. 

The suppresion of the reflection of the Alfven waves 
(the increase of ct with (La/)o ) is the main reason 
of the initial rise of the energy convcsion factor, ce(= 
LK,out/ {LaDq', Equation 19) shown in Fig. 3. A larger 
fraction of the Alfven waves transmits to upper regions, 
which directly leads to the larger final kinetic energy of 
the stellar wind. 

4-2.2. Atmospheric Structure 

We demonstrate simulated wind structures with differ- 
ent energy inputs. In Figure 5, we plot the time-averaged 
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Fig. 6. Comparison of the accelerations (forces per mass) 
by the magnetic pressure gradient (solid lines) and the gas 
pressure gradient (dashed lines) of an inactive case (Model 
Ba— 2a; left) and an active case (Model Cd-l-37; right). 

structures of the three cases tabulated in Table 2. The 
input Poynting energy of the Alfven waves from the pho- 
tosphere is scaled by (La/)o Br.ofoSvQ, with taking 
into accont the effect of the adiabatic expansion of the 
flux tubes. Then, from the flrst case (Model Ba— 2a) to 
the second case (Model CbOa) and the third case (Model 
Cd-l-37), the input energy flux increases by 8 times each. 

The top panel of Figure 5 indicates that the decrease 
of the density becomes slower with increasing the energy 
input from the surface. This is mainly because in the 
active cases the magnetic pressure dominantly supports 
the chromospheric density structure in addition to the gas 
pressure. Figure 6 compares the force (per mass) by the 
magnetic pressure gradient (solid lines) with the force by 
the gas pressure gradient (dashed lines) of the most in- 
active case (Model Ba— 2a; left panel) and the most ac- 
tive case (Model Cd-l-37; right panel) among the three 
cases. The left panel shows that in the inactive case the 
gas pressure gradient mostly dominates except around the 
transition region, 7' — rg ^ 10~^rQ. On the other hand, in 
the active case (right panel) the magnetic pressure gra- 
dient dominantly supports the density structure at and 
above the middle chromosphere, r > 2 x lO^^ro, because 
the input wave amplitude {Sv± oc 5B±) from the surface 
itself is larger. As a result, the cool chromospheric mate- 
rial extends to an upper region and the transition region 
is located at a higher altitude. The density of the above- 
located corona is also larger. 
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Fig. 7. Comparison of time-average Alfven wave properties 
of the three cases tabulated in Table 2. The line types are the 
same as in Figure 5. top: Energy flux of the Alfven waves. In 
order to remove the effectof the expansion of the flux tubes, we 
display FAr^//rQ = LA//47rrQ (see Equation 9 or 11), middle: 
Ratio of the outgoing and incoming (reflected) components of 
the Alfven waves, z'^/z'^. Bottom: Squread Elsasser variable. 

The change of the density structure directly affects the 
reflection of the Alfven waves. In the chromosphere the 
Alfven speed, va = Br/y/^np, generally increases with 
height because of the decrease of the density. The slower 
decrease of the density in the active case makes the 
change of va more gradual. Then, the reflection of the 
Alfven waves is not so severe in the active case. 

The wind velocity becomes slower with increasing the 
energy input because more mass is lifted up to the wind 
region, and then, the wind material cannot be effectively 
accelerated to higher velocities. The maximum temper- 
ature of the active case is highest, but the temperature 
drops faster because the thermal conduction is relatively 
ineffective in comparison with the adiabatic cooling in the 
high density circumstances 

Table 2 clearly indicates that the transmisivity (ct = 
(LA/)tc/ (LaDo) of the Alfven waves through the chromo- 
sphere increases with the input energy as shown in Figure 



Fig. 8. Kinetic energy, LK,outi (filled squares) and radiation 
loss, {L^f)tc (crosses) versus upgoing Poynting flux, (iA/)tc, 
at the top of the chromosphere of the cases with Brfifo = 1-25 
G. The diagonal line indicates the x = y relation and the 
horizontal line indicates the saturation level of -I/K,out ■ 

4. Accordingly, the kinetic energy luminosity, iK,out, 
and the mass loss rate, M, rapidly increase with (La/)o- 
The radiation loss, (LR/)tc, also rapidly increases with 
(La/)o: a-nd its dependence is faster than the dependence 
of ijK.out, which we discuss in §4.3 in terms of the satura- 
tion of the stellar winds. 

4.2.3. Wave Properties 

Figure 7 displays properties of the Alfven waves of the 
three cases shown in Figure 5 and Table 2. The top 
panel of Figure 7 compares the Alfven wave energy fluxes 
(Equation 9). Here we plot FAr^//r§ = LaJ l^irrl to ex- 
clude the effect of the adiabatic expansion in the super- 
radially open flux tubes. While the input Poynting flux 
((La/)o oc Bo/o^Vq) from the photosphere is different by 
8 times each, the differences of the wave energy fluxes 
in the atmospheres are more than one order of magni- 
tude. TFigure 12he sensitive dependence of the wave en- 
ergy flux in the atomosphere on the input wave energy is 
because of the suppression of the reflection of the outgo- 
ing Alfven waves as discussed so far. The middle panel 
of Figure 7, which compares z'^/z'^, the ratio of the re- 
flected waves to the outgoing waves, clearly exhibits that 
the wave reflection is suppressed for larger (iA/)o- 

The bottom panel of Figure 7 compares . This inter- 
estingly shows that the wave amplitudes of the outgoing 
component are comparable in these three cases although 
the energy fluxes are very different. The difference of Fa is 
mostly owing to the difference of the densities, in addition 
to the difference of the reflection efficiency. 

4.3. Saturation o/ix.out -Enhanced Radiative Loss- 

4.3.1. Energetics 

We examine the energetics of the stellar winds above the 
chromosphere by Equation (16). In Figure 8, we show the 
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Fig. 9. Energetics of the simulated stellar winds. From top 
to bottom, kinetic energy, LKouti radiation loss, [L^f )tc, 
and gravitational loss, {LQf)tc, are shown on upgoing 
Poynting flux, {L\f)tc, at the top of the chromosphere. The 
diagonal line in each panel indicates the x = y relation. In the 
panels for LK,out and (Lof )tc, the present solar level which is 
estimated by using M = 2 X 10~^* M^yr'^ and Vr = 400 km 
is also shown by the horizontal lines. In the top panel, 
the saturation level of Lk out fof each Br,ofo is indicated by 
the color arrows. 



wind kinetic energy luminosities, ^K.out (filled squares), 
and the radiation losses, {L^f )tc (crosses), with the net 
outgoing Poynting fluxes. {Lxf)tc, measured at the top of 
the chromosphere where T = 2 x 10"* K of the cases with 
Brflfa = 1-25 G. The fraction of the wind kinetic energy, 
iK,out/(.^A/)tc, is initially large, with roughly half of 
{Lxf)tc transferred to the stellar wind. LK.out increases 
almost linearly with {LA,f)tc, but eventually saturates for 
large [Lxf)tc- Some cases even show considerably smaller 
-^K.out than other cases with similar {L p^f )tc- The satu- 
ration of the wind kinetic energy is a consequence of the 
enhancement of the radiation loss. {L^f )tc increases with 
{Lxf)tc much faster than the linear dependence, because 
the wind density increases as a result of strong driving 
as we increase the input energy. In the optically thin 
plasma, the radiation loss is proportional to p^, and then, 
(LR,/)tc rapidly increases with (LA/)tc. The fraction of 
the radiation loss, (iR/)tc/(iA/)tc, which starts from a 
much smaller value at small (LA/)tc than iK,out/(^A/)tc, 
is approching to sa unity^ ((iR/)tc = (LA/)tc) for large 
{Lxf)tc- In this saturated state, the Poynting flux passing 
through the top of the chromosphere is mostly radiated 
away in the upper regions. As a result, the fractions trans- 
fered to the wind kinetic energies are quite small in these 
cases, which is observed as the saturation or drop of -/jK,out 
with increasing (-/jA/)tc. 

Figure 9 displays each term of Equation (16), the ki- 
netic energy luminosity, iK.out (top panel), the radia- 
tion loss, {Ljif)tc (middle panel), and the gravitational 
loss, (LG/)tc (bottom panel), of the stellar winds with 
the net outgoing Poynting flux, {Lxf)tc, measured at the 
top of the chromosphere. Figure 9 clearly shows that 
the trend -the saturation of ^K.out as a consequence of 
the enhanced {L^f)tc - is universal for the other values 
of Brfifo- Furthermore, the saturation level of the ki- 
netic energy luminosity, Lk, out, sat, is positively correlated 
with Br^ofo as shown in the top panel, where Lk, out, sat 
is taken from the case which gives the largest LK,out for 
each Brfifo (Table 3). The nonlinear dissipation of the 
Alfven waves can explain the dependence of the satura- 
tion level on B^.o/o- Since Br =i?r,o/o''o/^^i Br.ofo deter- 
mines the Alfven speed, va, in the wind acceleration re- 
gion (1.5ro ^ r ^ lOro)). The nonlinearity of Alfven waves, 
dv±/vA, is an important factor which controls the dissi- 
pation of Alfven waves, and larger 5v±/va leads to faster 
dissipation (e.g. HoUweg 1973; Suzuki & Inutsuka 2006). 
Then, in cases with smaller Br.ofo, the nonlinearity be- 
comes relatively larger because of the smaller va in the 
wind acceleration regions (Suzuki 2012). In these cases, 

^ In some cases, the radiative losses exceed the injected Poynting 
fluxes from the top of the chromosphere, (I'R/)tc/(i'A/)tc > 1, 
which seems to break the energy conservation relation. In ad- 
dition to the Poynting fluxes, however, there are contributions 
from perturbations of gas pressure (Eq.6) to the energy injec- 
tions, which we do not take into account in Equation (16). The 
gas pressure perturbations are due to sound waves nonlinearly 
generated from Alfven waves (Kudoh & Shibata 1999; Suzuki 
& Inutsuka 2006). Although the energy inputs from such sound 
waves are smaller than the energy inputs from the Alfven waves, 
they are not negligible in active cases. 
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Model 


Svo 


Br,o(kG)x/u 


hi (ro) 


(LAf)o 


(iA/)tc 


CT 


M (A/oyr-i) 




^K.out 


(iR/)tc 


Ac+27 


2.68 


2/6400 


0.1 


1.2 X lO^'-' 


1.7 X 10^» 


0.14 


1.1 X 10"^'' 


226 


1.8 X 10^' 


2.4 X 10^' 


Bc+4/3 


5.35 


2/3200 


0.03 


9.8 X lO^'-* 


1.1 X lO^'-* 


0.11 


5.6 X 10"^-' 


278 


1.4 X 10^» 


1.8 X lO^'-* 


Cd+37 


3.79 


4/3200 


0.1 


9.8 X 10^" 


1.6 X lO^'* 


0.17 


8.5 X 10-^'^ 


313 


2.6 X 10^» 


2.2 X 10^» 


Df+37 


3.79 


16/6400 


0.1 


2.0 X 10^" 


5.7 X 10^" 


0.29 


2.5 X 10"^- 


373 


1.1 X lO^'-* 


7.2 X 10^** 


Ef+47 


5.35 


16/3200 


0.1 


7.8 X 10'^" 


2.1 X lO'^" 


0.26 


1.0 X 10""^ 


351 


3.9 X 10^^ 


4.5 X lO^'-* 


Ff+57 


7.57 


16/1600 


0.1 


3.1 X lO'^^ 


8.2 X 10-^" 


0.26 


4.3 X 10"^^ 


317 


1.4 X 10'^" 


2.9 X 10-^" 



Table 3. Saturated cases for Bj.fi fo = 0.3125, 0.625, 1.25, 2.5, 5, 10(G) from top to bottom, shown in Figure 10. The unit of &vq and 
fr.out is km s~^, and the unit of L is erg . cx is the surviving fraction of the Alfvcn waves from the photosphere to the top of 
the chromosphere (Equation 15). 




Br,ofo (G) 

Fig. 10. Upper panel: Kinetic energy, Lk, out. sat, of the sat- 
urated cases for the respective B^-fifo with i?r, o/o- The hne 
is the best fit function. Equation (21). Lower panel: Mass 
loss rate, Mgat of the saturated cases. The line is the best fit 
function. Equation (22). 

the Alfven waves dissipate faster at lower altitudes where 
the densities are higher. As a result, more wave energy 
is transfered to the radiation losses rather than the wind 
kinetic energy. On the other hand, in cases with larger 
Br.ofo the Alfven waves can travel to higher altitudes and 
effectively drive the stellar winds, which gives the higher 
saturation levels of LK,out- 

We show the wind kinetic energy, Lk, out, sat, and the 
corresponding mass loss rate, Afgat , of the saturated cases 
(Table 3) in Figure 10. Both Xk. out. sat and Mgat are well 
fitted by power-law functions (solid lines in Figure 10), 



iK,out,sat = 2.05x 1028crgs-i(S,^o/o)''^(^^^ ,(21) 
and 

Ms.t = 7MxlO-'^Meyv-\Br^ofo)'-'' '(22) 

where we take into account the tq and M dependences 
in order to extend to general solar-type stars and directly 
compare with Equations (19) and (20). Both relations 
give the similar dependences on ( i?r, 0/0) • 

The comparison between Equations (19) and (21) gives 

CE(H'>l-toc(i?,,o/o)°-'' (23) 

for the saturated cases, where we neglect the weak depen- 
dence on the photospheric density, po- In order to under- 
stand this relation, let us start from the dependence of ce 
on Svq in a fixed (i3r,o/o)- In the unsaturated state with 
small {5vq), an increase of (Svq) enhances ce, because the 
reflection of the Alfven waves is suppressed (§4.2). When 
one further increases (Sv^) to the saturated state, ce starts 
to decrease on account of the radiation loss as discussed 
in this section and finally cb{Svq) itself decreases. At the 
saturation, ce ~ 0.05, which is nearly independent from 
{Br,ofo) (Figure 3). Equation (23) indicates that {5vq) 
that gives the saturation of the stellar wind is roughly pro- 
portional to (i?r,o/o) (we assume 0.84 « 1.). In cases with 
large (i?r,o/o), large amplitude Alfven waves can propa- 
gate to higher locations owing to larger wa and accordingly 
the smaller nonlinearity of the Alfven waves as discussed 
so far. 

Equations (21) and (22) are the maximum kinetic en- 
ergy and mass loss rate for given (i?j.,o/o)- The stan- 
dard case for the current solar wind (Model CbOa; 
Table 2) gives only 8% of the maximum kinetic energy 
(-£'K,out/iK,out,sat = 0.08). Therefore, by increasing Svq at 
the photosphere, the stellar wind kinetic energy could be 
raised up to 12.5 times but cannot exceed it. In order to 
get further large kinetic energy, a fiux tube with larger 
{Br.ofo) is required. 

An important point regarding Equations (21) and (22) 
is that the saturation level is not determined by only Brfl 
but Brflfo; open fiux tubes need to occupy a sufficient 
fration of the stellar surface {i.e., large /o) in addition to 
large field strength at the photosphere. This is also con- 
sistent with the explanation of the saturated stellar winds 
from the change of the magnetic topology that closed mag- 
netic loop structure, which is expected to dominanty cover 
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Model 


Svo 


Br,o(kG)x/o 


hi (ro) 


{LaDo 


(iA/)tc 


CT 


M (Moyr-i) 


t'r,out 


LK.out 


(iR/)tc 


Dd+37 


3.79 


4/1600 


0.1 


2.0 X lO-'" 


1.2 X 10^" 


0.061 


3.3 X 10-^'' 


733 


5.6 X 10^» 


4.9 X 10^' 


Dc+47 


5.35 


8/3200 


0.1 


3.9 X 10-'" 


1.1 X 10-'" 


0.28 


6.0 X 10"^^ 


145 


4.0 X 10^» 


4.7 X lO^'-* 


Df+57 


7.57 


16/6400 


0.1 


7.8 X 10^" 


1.6 X 10-^" 


0.21 


4.0 X 10-^^ 


191 


4.6 X 10^» 


1.5 X 10^" 



Table 4. Input parameters and outputs for the cases shown in Figure 11. The unit of i5?;o and Dr,out is km s ^, and the unit of L 
is erg s~^. Ct is the surviving fraction of the Alfven waves from the photosphere to the top of the chromosphere (Equation 15). 
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Fig. 11. Comparison of the time-averaged wind structures 
of three cases tabulated in Table 4. The less active case 
(Model Dd-l-37), the intermediate case (Model De-l-47), and 
the active case (Model Df-l-57) are shown in dotted, solid, 
and dashed lines. From top to bottom, the densities, radial 
velocities, and tcmcpraturcs arc compared. 

the surface of active stars, inhibits the stellar winds (Wood 
et al. 2005). A large fraction of closed loops corresponds 
to a small /o, and then, stellar winds arising from such 
conditions will also become weak in our simulations. 

4-3.2. Atmospheric Structure 

We further study the saturated state of the stellar winds 
by examining some simulated wind structures. Figure 11 
compares the three cases, summarized in Table 4. We 
select runs with the same Sr.o/o = 2.5 G, but different 
Svq. From the first case (Model Dd+37) through third 



case (Model Df+57), the input energy flux (oc BrfifoSv^) 
from the photosphere increases by twice each. 

From the first case (Model Dd+37) to the second case 
(Model Dc+47) , the atmospheric structure drastically 
changes. The density of the second case is ^10 times 
larger in the chromosphere and corona and '--^lOO times 
larger in the stellar wind region. On the contraty, the wind 
speed becomes very slow because the dense wind cannot 
be effectively accelerated. As a result of the high density, 
the radiative loss {{L-nf )tc in Table. 4) of the second case 
is nearly 100 times larger than that of the first case, even 
though the input wave energy is only twice. Although 
the mass loss rate (Af c>c pvr in Table. 4) is also larger 
due to the enhanced density, the kinetic energy luminosity 
(^K.out oc pwj? in Table. 4) slightly decreases because of 
the low velocity, which is recognized as the saturation of 
-f'KjOut in the top panel of Figure 9. While the maximum 
temperature of the second case is higher than that of the 
first case, the temperature decreases more rapidly because 
the thermal conduction is relatively ineffective against the 
adiabatic cooling in the denser corona. 

The comparison between the second (Model De+47) 
and third (Model Df+57) cases exhibits an example of a 
saturated state. The overall density and velocity struc- 
tures are quite similar. Then, both M and iK.out from 
these two cases are similar; the mass loss rate also satu- 
rates at this stage. By a close look one can see that the 
third case gives slighly larger density from the chromo- 
sphere to the low coronal region, which leads to the ^ 3 
times larger radiative loss as shown in Table. 4. The max- 
imum temperature of the third case is smaller than that 
of the second case because of the enhanced radiative loss. 

5. Discussions 

5.1. Comparison with Wood et al. 

In this section, we compare our results with the obser- 
vation by Wood et al. (2002; 2005). As a key presentation 
in their papers, they plotted the mass loss rates with the 
X-ray fluxes of the observed solar- type stars; we here put 
our results in the same diagram. The mass loss rates ob- 
tained in our simulations can be directly compared with 
the observations. On the other hand, the X-ray fluxes 
of our simulations considerably underestimate the actual 
X-ray fluxes because our simulations only treat open flux 
tubes although it is expected that the observed X-rays are 
mainly from closed loops. From an observational point of 
view, the UV flux from solar-type stars is well-correlated 
with the X-ray flux (Ayres 1997). Based on these con- 
siderations, we use the radiation flux from the gas with 
T > 2 X 10^ K. 



12 



Suzuki et al. 



B,ofo(G)=10 □ 
5 o 
2.5 
1.25 



0.625 

0.3125 □ 

Wood+ 2005 (MS) • 

Wood+ 2005 (EV) X 




2(LRt),/47iro or Fx obs(erg cm" s; 



V □ 

o 



P 5 



X 



10" 



10' 



10" 



10" 



10' 



10° 



2(LRf),/47ir^ or Fx obs(erg cm"^s) 



Saturation Level ^ ^ " 

for each jfp 5 ^ 

□ O • V 



1.25 * 
0.625 - □« 



0.312 5 , ^ 



o 0° 00 



X 



-X- 



10" 



10' 



10" 



10° 10' 

2(LRt),/47ir^ or Fx obs(ei'g cm"^s) 



10° 



Fig. 12. Comparison of the simulation results (color sym- 
bols) with the observed data by Wood et al.{2002; 2005; 
big black symbols; circles indicate main sequence stars and 
crosses are evolved stars.) The horizontal axis displays the 
X-ray flux (for Wood et al. data) or the radiation from open 
flux tube regions with T > 2 X 10" K multiplied by Cr = 2 
(for our results; see text for detail). From top to bottom we 
compare the mass loss rates (oc pwr.out), the ram pressures 
(oc P'^^out)' ^'^'^ kinetic energy luminosities (oc P'f^out) 
divided by the stellar surface areas and normalized by the 
present-day solar level (M = 2 X lO~l*M0yr~l and v = 400 
km s~^). Then, the observational data points are plotted at 
the same positions in the three panels, because Wood et al 
(2002; 2005) adopt the constant v =400 km s'^ for their hy- 
drodynamical simulations to estimate the observational M. 
The color arrows in the bottom panel indicate the saturation 
levels of I/K,out for the respective i?r,o/o- 



In order to compare our results with the observed X- 
ray fluxes, Fx, we need to estimate a conversion factor, 
Cr, from the radiation flux from the gas with T > 2 x lO'' 
K in the simulated open flux tubes to the actual X-ray 
flux averaged over the entire stellar surface, where Cr is 
explicitly in an equation, 

(iR/)tc 



Fx 



Here, 



47rrg 



(24) 



, ; is the radiation flux from the gas with 

r > 2 X 10* K in the open flux tubes. While it does not 
include the X-ray from closed loops as discussed above, 
it includes the radiation in the UV range in addition to 
the X -ray; Cr could be either larger or smaller than unity. 
Because we do not have a quantitatively accurate way to 
estimate Cr, we simply choose a value of Cr to give a rea- 
sonable fit to the observation. In Figure 12, we show our 
results with Cr = 2 in comparison with the observed data. 
From top to bottom, we compare the mass loss rates, M, 
ram pressures, Mv^oxa, £md kinetic energy luminosities, 

M , divided by the surface areas of the stars and 
normalized by the current solar level. The three panels 
essentially compare pwr.out, P^^r.out: ^-^d pfj^ o^f We use 
M© = 2 X 10"^^Mq yr"i and f© = 400 km s"^ as the nor- 
malizations in Figure 12 according to Wood et al. (2002; 
2005) in which they estimated M by comparing the ob- 
served Lya absorptions with their hydrodynamical sim- 
ulations of the astrosphercs with a constant wind speed 
= 400 km s~^ for the different stars. Since they assume 
the constant wind speed, we put the observed data points 
at the same positions in the three panels. 

One can recognize that our simulations exhibit more 
clear saturations in pv^ and pv^ . On the other hand, 
pVr.avA. shows an increasing trend rather than a satura- 
tion. This is mainly because active cases with large ra- 
diation fluxes give systematically slower wind velocities 
(Fig. 5). The Lya absorptions in the hydrogen walls of 
astrospheres are correlated with the ram pressures of stel- 
lar winds (Wood & Linsky 1998). In this sense, the ram 
pressure (the middle panel) is a more physically plausi- 
ble variable. In this panel the simulation results explains 
the observed distribution quite well except for the evolved 
stars. 

In realistic situations, the average Sr.o/o of the open 
flux tubes of s sun-like star would change with time. 
Comparison of the observed values with the simulation re- 
sults in the middle or bottom panels of Figure 12 enables 
us to raise a following speculative scenario on the time 
evolution of the stellar wind (Figure 13). Stars evolve in 
general from the right to the left in Figure 12 . The most 
active star among the observed data at Fx = 1.3 x lO^erg 
cm~^s~"'^ is well explained by a small value of i?r, 0/0 < 1 
G. Probably, this star is mostly covered by closed loops 
and the filling factor of open flux tubes is very small, say 
/o ^ 1/ 10000. Then, even if the photospheric field strength 
is strong, e.g. 5^,0 ~ 10 kG, i?r,o/o becomes small < 1 G, 
which gives the low saturation level of the stellar wind. 
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Fx ~ 10' erg cm 
Brfi 10 kG,fo<j^ 

Saturated 

M < lOMo 



AAAAA 



AAAAA 

Fx ~ lO^erg cm ^ ^ 10^ erg cm ^ 



Brfi ~ 5 kG, /o ~ 
Saturated 
M > lOOM© 



1000 



Brfi - 1 kG, /o - jooo 
Unsaturated 



Fig. 13. Speculative scenario of the time evolution of properties of open flux tubes based on the observation by Wood et al. (2002; 
2005). (see text) 



With the stellar evolution, the magnetic topology is re- 
laxed so that /o increases to ^ 1/1000. At this stage with 
Fx ^ lO^erg cm~^s~^ , the stellar wind flux increases owing 
to the increase of /o, although stars are still in the satu- 
rated state. In the subsequent evolution, Brfi decreases 
with time and stars are deviated from the saturated state. 
The stellar wind kinetic energies eventually settle down to 
the current solar level. 

5.2. Faint Young Sun Paradox 

The luminosity of sun-like stars gradually increases dur- 
ing the main sequence phase. Based on the standard stel- 
lar evolution calculation, the luminosity of the sun is 20- 
30% smaller at early times (Gough 1981). As a result, 
the surface temperature on the earth is expected to be be- 
low the freezing temperature before 2 billion years ago. 
On the other hand, liquid water already existed at that 
time on the earth as well as on the mars (Sagan & Mullen 
1972; Karhu & Epstein 1986; Feulner 2011), which looks 
contradiction, known as the 'faint young Sun paradox'. A 
high level of green house effects (Owen et al. 1979; Kasting 
1997); and modification of the solar model (Whitmire et 
al. 1995; Minton & Malhotra 2007) are the two major 
possibilities which can solve the problem. We here dis- 
cuss effects of the strong mass loss of the Sun (Sackmann 
& Boothroyd 2003). 

As already discussed in §4.3, the saturation level of 
the stellar winds is determined by Sr,o/o- The top panel 
of Figure 9 and Figure 12 indicate that stars with large 
Br,ofo (~ 5-lOG) could drive ~ 500 — 1000 times stronger 
mass loss with the X-ray flux in the observed range (Fx < 
2 X 10^ erg cm~^s~^), although such a large mass loss rate 
has not been observed so far. If such a strong mass loss 
actually continues ^ 1 billion years, it contributes to the 
decrease of stellar masses. 



1000 times of M during 1 billion years gives the mass 
loss of O.O2M0. If this is true in the Sun, the initial mass 
is 2 % larger than the present Sun. An increase of the stel- 
lar mass influences the radiation flux on planets around 
a star in two ways. First, the stellar radiation itself is 
larger; in sun-like stars, the luminosity, L, can be scaled 
by Loc AF'-'^^ (Kippenhahn & Weigert 1990; Whitmire et 
al. 1995). A 2% larger M gives w 10 % larger L. Second, 
the orbital radius, Op, of a planet shrinks in a manner 
to keep A/cp constant (Whitmire et al. 1995; Minton & 
Malhotra 2007). Then, the radiation flux (cx 1/a^) on 
the planet increases by w 4 %. As the sum of these two 
effects, the radiation flux increases « 15 %, which could 
considerably compensate the early faint sun and avoid the 
freezing temperature on the earth. 

In this paper, we are focusing on the quasi-steady com- 
ponent of the stellar winds from open flux tubes during 
the main sequence phase. In active sun-like stars, coronal 
mass ejections (CMEs) , which probably involves magnetic 
reconnections of closed loops (e.g. Ohyama & Shibata), 
might significantly contribute to the mass losses. Aarnio 
et al. (2012) estimated mass loss rates by CMEs from 
pre-main sequence stars by extrapolating the observed re- 
lation between X-ray flare energy and CME mass in the 
present Sun. Although their main target is not CMEs 
during early main sequence phases, which is crucial in the 
faint young Sun paradox, we can infer by a simple inter- 
polation that the relative contribution of CMEs to the 
total mass loss is larger at earlier times than the present 
contribution (~ 10%). 

5.3. Extended Variable Chromosphere 

We have discussed in §4.2 that, when increasing the 
input wave energy from the current solar level, the chro- 
mospheric materials are lifted by the magnetic pressure 
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Fig. 14. Comparison of snap-shot wind structures at 
t = 9.296 X 10^ (s) (solid) and t = 9.707 x 10^ (s) (dashed) 
for Model Dd+47. The dotted lines are the time-average be- 
tween t = 3.328 X 10^ (s) and 1.331 X 10*^ (s). 



of the Alfven waves (§4.2.2 and Figure 6). An extended 
chromosphere is a universal feature in our simulations 
with relatively large energy inputs. Another important 
aspect is that the atmospheres behave more dynamically. 
In particular, the transition regions move up and down 
because of the thermally unstable region (T ^ 10^ K) 
of the radiative cooling function (Landini & Monsignori- 
Fossi 1990; Sutherland & Dopita 1993). Figure 14 shows 
snapshot structures (solid and dashed lines) of an active 
case (Model Dd+47) with Svq = 5.35 km s~\ Br^ = 4 
kG, fo = 1/1600; and h\ = O.lro, in comparison with the 
time-averaged wind structure (dotted lines). In the time- 
avcragcd structure the temperature gradually increases 
from T = IQ-* K at r w O.Olro to T = 10^ K at r w 0.15ro. 
However, one can recognize from the snapshots that this 
is simply a result of the long-time average. The snapshot 
structures at the two different times show that the sharp 
transition region moves up and down; at t = 9.296 x 10^ s 
(solid lines) the chromosphere is extended up to r w 0.2ro, 
white at t = 9.705 x 10^ s (dashed lines) the corona is get- 
ting down to r « O.Olro with a part of dense coronal gas 
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Fig. 15. Location of the top of the chromosphere, rtc, on 
radiation flux of each run. The filled symbols distributed 
in the upper half of the panel indicate the maximum heights 
during the simulations, and the open symbols near the bottom 
of the panel indicate the minimum heights. In some cases 
with large radiation flux, (LR/)tc/47rrQ, the minimum of rtc 
is below the displayed region in the panel. 

falhng down (Pinto et al. 2009). 

An extended chromosphere is actually observed. 
Recently, Czesla et al. (2012) observed the atmosphere of 
a planet-hosting young star, CoRoT-2a, by using the tech- 
nic of the Rossiter-McLaughrin effect. They reported that 
it has a very thick chromosphere with height of 10-20% of 
the radius, which is much larger than the chromosphere 
of the present Sun which is only up to < 1% of the ra- 
dius. Moreover, they observed an asymmetry in the chro- 
mospheric Rossiter-McLaughrin effect, which might arise 
from inhomogeneities of the chromospheric structure. 

Figure 15 shows the height of the chromosphere of each 
simulation run. In the figure, we plot both the maximum 
(filled symbols) and minimum (open symbols) values of 
rtc at T = 2 X 10"* K during the time evolution of each 
simulation. In cases with small radiation flux, the dif- 
ferences between the maximum and minimum heights are 
small and rtc is typically O.Olro, which is consistent with 
the present-day Sun (e.g. Imada et al.2011). As the ra- 
diation flux increases, the maximum height increases and 
the minimum height decreases, namely in active cases the 
transition regions largely move up and down. The mini- 
mum heights in cases with relatively large radiation flux 
(in the right part of the figure) are considerably small, 
''tc — ''o < 10~''ro(= 70 km). Moreover in some cases the 
minimum heights are in the outside of the displayed range. 
In these cases, the temperature increases to T > 2 x IC* K 
from the next mesh (sa 10~^ro = 7 km) to the photosphere 
at some instants. 

5.4- Stellar Rotation 

It is expected that stellar rotation is an impor- 
tant ingredient in determining stellar magnetic activities 
(Skumanich 1972; Ayrcs 1997; Giidel et al. 1997; Linsky 
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et al. 2012). Probably related the steep differential rota- 
tion owing to the fast rotation of young active stars (Hotta 
& Yokoyama 2011), the magnetic field is effectively ampli- 
fied in fast rotating solar-type stars (Gondoin 2012). We 
take into account the effects of the strong magnetic fields 
on the stellar winds by incorporating the wide ranges of 
the parameters on the magnetic field. 

In addition to the generated magnetic field strength, 
the stellar rotation directly affects the dynamics of stel- 
lar winds through the magnetocentrifugal force (Weber & 
Davis 1967; Cohen et al. 2010). The angular momentum 
of a star is outwardly transported by the magnetorotating 
wind, which is the main cause of the deceleration of the 
stellar rotation (Matt et al. 2012). We do not take into 
account this effect in our simulations. We here briefly dis- 
cuss possible roles of the magnetocentrifugal force in the 
dynamics of the stellar winds. 

In general, stronger magnetic field and faster rotation 
affects the dynamics of stellar winds. Magnetic energy per 
mass in magnetorotating winds (Weber & Davis 1967; see 
also §9 of Lamers & Cassinelli 1999) is written as 

= -. ~ — ; (,25j 

where 51 is the rotation frequency of a star, is the 
azimuthal component of the background magnetic field 
{BrB^ < because of winding field lines), and we here 
consider the problem in the outer region where the super- 
radial expansion of the flux tubes already finishes, f — 1. 
For the second transformation in Equation (25), we have 
used the relation for the geometry of magnetic field lines 
under the steady-state condition, 

B^ _ v^-rn ^ rn ^^g^ 

where is the azimuthal component of velocity. If eb is 
larger than kinetic energy per mass, the magnetocentrifu- 
gal force dominantly works and the wind velocity will be 
faster than obtained in our simulations. Then, we can use 
a nondimensional parameter. 



47rr2pw4/2 



iBrM)'n^ 

Mvf/2 



(27) 



which measures the importance of magnetocentrifugal 
force^. The effect of magnetocentrifugal force becomes 
important in proportion to magnetic energy and rotation 
energy but in inversely proportion to a mass loss rate. 

We can estimate Csn from a typical case of our simu- 
lations for active stars. 



Cen - 0.57 



5G 



M 



—) (-) 

MqJ yReJ l^lOOMg 



400km s-i 



O.IOk 



,(28) 



^ What wc arc doing here is essentially the comparison between 
the Michel velocity, vj^, (Michel 1969) and the escape veloc- 
ity, Vcsc,o- Neglecting numerical factors with an order of unity, 
CsQ = 1 corresponds to vm ~ Vescfi- 



where M© = 2 x 10 ^^AlQ-yr ^ is the mass loss rate of 

the current Sun, and ~ \ is the Kepler (breakup) 

V ''o 

velocity at the stellar surface. Equation (28) is normal- 
ized by 10 % of 51k, which gives the dependences on M 
and rp. The rotation of the current Sun is VLq « 0.005flK- 
Then, Equation (27) indicates that the effect of the mag- 
netocentrifugal force is small if the rotation frequency of a 
star which gives 100 times of M is smaller than 20 times of 
the current solar rotation frequency. In stars with rotation 
period ^ 1 day the magnetocentrifugal acceleration is im- 
portant and the termical velocity becomes faster, whereas 
it also depends on the mass loss rate. 

6. Summary 

We have investigated how the properties of the stellar 
winds from sun-like stars are determined from the mag- 
netic fields and velocity perturbations at the photosphere. 
We performed the 163 models of the MHD simulations in 
the range of ~ 4 orders of magnitude of the Alfven wave 
energy fiux from the photosphere. We examined the prop- 
erties of the stellar winds of these simulations mainly from 
an energetics point of view. 0.1-10%, with the typical 
value, ^ 1%, of the input Poynting flux of the Alfven waves 
from the photosphere is finally transported to the kinetic 
energy of the stellar winds. We derived the scaling rela- 
tions that dermine the kinetic energy luminosity and mass 
loss rate of the stellar winds from the properties at the 
stellar photosphere (Equation 19 & 20). These relations 
can be used to estimate the mass losses from solar-type 
stars (Holzwarth & Jardine 2007; Cranmer & Saar 2011). 

We can classify the simulated stellar winds into the un- 
saturated regime and the saturated regime. When the 
input wave energies are small, the stellar winds arc in the 
unsaturated state; an increase of the input wave energy 
directly leads to strong stellar winds. In this regime, the 
refiection of Alfven waves dominantly controls the ener- 
getics and dynamics of the stellar winds. In the chromo- 
sphere and low corona, the density rapidly decreases with 
height. Consequently, the Alfven speed rapidly increases, 
which makes Alfven waves effectively refiected because 
of the deformation of the wave shapes. The reflection 
is more severe for a smaller wave energy input because 
the density decreases more rapidly without the support 
from the magnetic pressure (oc 5B\ /Stt) associated with 
the Alfven waves. In extreme cases, more than 99% of the 
input Poynting energies are reflected back before reaching 
the top of the chromosphere. The mass loss rates of these 
cases are as small as ~ 1% of that of the current sun. 

With increasing the input wave energy, the density de- 
creases more slowly in the chromosphere owing to the sup- 
port by the magnetic pressure of the Alfven waves. The 
chromosphere of the active cases is time-dependently ex- 
tended to 10 - 20 % of the radius, and the transition region 
moves up and down (§5.3) mainly because of the ther- 
mally unstable region {T ^ 10'"' K) of the radiative cooling 
function of the gas with the solar abundance (Sutherland 
& Dopita 1993). According to the slow decrease of the 
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density in the chromospheres, the Alfven speed changes 
more gradually, which suppresses the reflection of the 
Alfven waves. When one increases the input wave energy, 
the kinetic energy and the mass loss rate of the stellar 
winds increase more rapidly than the linear dependence 
on the input energy. In cases with the large energy inputs 
for active stars, more than 10% of the input wave energy 
from the photosphere reach the top of the chromosphere. 

For the large wave energy inputs, the stellar winds are 
in the saturated state; an increase of the input energy does 
not raise the kinetic energy of the winds but enhances the 
radiative loss. Increasing the input wave energy from the 
unsaturated state to the saturated state, the densities of 
the transition region and corona monotonically increase, 
which gives larger kinetic energy and mass loss of the stel- 
lar winds. The increase of the densities also raise the ra- 
diative cooling. Since the radiative loss is proportional 
to in the optically thin plasma, the increase of the ra- 
diative loss is faster than the increase of the wind kinetic 
energy with increasing the input wave energy from the 
photosphere. In the saturated state, most of the Poynting 
energy at the top of the chromospheres is transferred to 
the radition in the transition regions and coronae, and 
a tiny fraction of the energy is available for the kinetic 
energy of the stellar winds. 

The saturation level has the power-law dependence 
on Br.ofo: -^K, out, sat CK {Brfifof -^^ (Equation 21), which 
can be understood from the nonlinear dissipation of 
the Alfven waves. Brflfa determines the magnetic field 
strength in the upper corona and wind region where the 
super-radial expansion of the flux tubes already finishes. 
For small i?r,o/oi the nonlinearity of the Alfven waves, 
SBj^/Br becomes systematically larger, which results in 
faster dissipation. Then, more wave energy dissipates at 
lower altitudes where the density is higher, and mostly es- 
capes by the radiation loss. For large Br^fo, the opposite 
explanation can be done; more wave energy remains up 
to higher altitudes and contributes to driving the strong 
stellar winds. 

The positive correlation of the saturation level with 
Br.ofo indicates that in order to drive strong stellar winds 
the filling factor of open flux tubes should not be so small, 
as well as the surface magnetic field should be sufficiently 
strong. This is consistent with the explanation that closed 
magnetic loops inhibit stellar winds from active stars. 
Based on our results, we introduce a speculative scenario 
for the evolution of the stellar winds from sun-like stars 
(§5.1; Figure 13). At very early time, the surface of a star 
is mostly covered with closed magnetic field (small /o), 
and then, the mass loss rate is not so large because of the 
low saturation level, even though the surface field strength 
{Br,o) is large. The magnetic field structure is eventually 
relaxed so that open flux tubes occupy a larger fraction 
(/o increases) and the mass loss rate increases. After that, 
the mass loss rate gradually decreases to the current solar 
level with decreasing magnetic field strength. 
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